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Abstract 

This article gives a brief overview of a "fundamental solution" based energy-norm harmonic inter- 
polation approach for two half-space settings of interest: the upper-half R n plane, where fundamental 
solutions satisfy Laplace's equation, and the upper-half complex plane, where simple poles can be re- 
garded as the fundamental kernels of interest (versus fundamental solutions). It is also pointed out that 
the approach can handle higher-order pole fits, as well as logarithmic source fits, in the complex setting 
and that it can handle higher-order multipole fits in the general real R n setting. Higher-order multipoles 
in the real R™ half-space setting are of particular interest since fits based on a commonly used type of 
radial-basis functions, which are known as inverse multiquadrics, can be reinterpreted as multipole based 
interpolations with the formalism presented here; moreover, this reinterpretation has certain significant 
previously unrecognized theoretical implications. 

Although this article can be read independently from the articles that preceded it in the series, 
a brief look at the content of these preceeding articles is useful here. In the first article a simple 
interpolation approach that satisfies a minimum energy-norm condition was presented for representing 
analytic functions in the exterior of a complex unit disk by a linear combination of simple poles located in 
the disk's interior. Corresponding interpolation techniques for the exterior of real unit disks and for the 
interior of real and complex unit disks were also presented. Relationships of this approach to Bergman and 
Szego kernel interpolation theory were also addressed. The second article gave a corresponding technique 
for using reciprocal-distance fundamental solution basis functions to perform harmonic interpolation over 
the interior of an R 3 sphere. This R 3 technique is based on closed-form weighted Dirichlet integral inner 
products, which means that when source positions are specified in the exterior of the harmonic sphere 
a closeted-form linear equation set results for the source parameters and that the resulting solution 
minimizes the weighted field energy of the modeling errors. It was also pointed out in this second article 
that a similar formalism exists for the exterior of sphere and for R 3 half-space, where the Dirichlet 
integral itself is directly minimized. The primary associated functional analysis setting for all of these 
developments is labeled Dirichlet integral dual-access collocation-kernel space (DIDACKS). Both previous 
articles showed that the basic DIDACKS approach can be easily extended to handle higher-order (multi-) 
poles. Finally, this second article also presented a proof that not only is the (weighted) Dirichlet integral 
norm of the modeling difference minimized, but also that of all candidate interpolating functions, the 
DIDACKS solution is the one with a minimum norm. 

This article derives analogous results for complex half-space and then generalizes the R 3 half-space 
results to R n half-space. This means that closed-from expressions are presented that yield interpolating 
fits, which directly minimize Dirichlet integral expressions in either C and R 3 half-space. Moreover, the 
DIDACKS solution satisfies a second minimization condition, which states that any interpolating function 
that matches the prescribed point data values must have a field energy content that is at least as large 
as that of the DIDACKS solution. This also means that when the above reinterpretation of inverse 
multiquadric radial-basis function fits as multipole DIDACKS fits can be carried out, these radial-basis 
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function fits also simultaneously satisfy two energy-norm minimization conditions, which in itself implies 
a certain level of numerical robustness. 

Finally, this article also shows that the same sort of development can be carried out for R n half-space 
surface integral norms where the square integral of the modeling error itself is directly minimized over 
the half-space boundary (hyper) plane. This development allows for a reinterpretation of radial basis 
fits similar to that just described; moreover, it also yields a mathematically motivated approach to the 
problem of half-space downward continuation. 

Key words: complex half-space, interpolation, Dirichlet form, radial basis function, Laplace's equation, 
complex poles , inverse problem, Dirichlet form, point sources, multipole, potential theory, 
fundamental solutions 
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1 Introduction 

This article develops a simple closed-form interpolation approach that minimizes the associated field energy 
over the half-space region of interest by using a linear combination of poles in C, or fundamental solution 
basis functions in R™, as interpolating functions. This means that in both settings, the resulting formalism 
can simultaneously be considered as an interpolation technique and as an approximation technique. 

First consider the approximation aspects of the resulting formalism. In the complex setting, where 
z = x + iy, the formalism approximates a given function f(z) that is analytic in the upper-half plane 
H c := {z E C | y > 0} by a finite set of spatially distributed simple poles in the lower-half plane H ± := {z E 
C | y < 0} so that the specific approximating (and interpolating) form becomes 

k— 1 

with \ik E C and Zk E . 

Likewise, in the real R" setting, where X = (aq, X2, %3, x n ) T , if f(X) is harmonic in the Tri- 
dimensional upper-half plane H" := {X E R™ | x n > 0}, then it is to be approximated by a finite set of spa- 
tially distributed fundamental solutions that have sources in the lower half-space H" := {X E R" | x n < 0}, 
so that the approximating (and interpolating) form becomes 

y^^JM^ (2) 

with fi k E R and I[ e H". (Other various notational conventions were explained in |H|S].) 

For a given /, when the values of the source locations are specified, the usual approximation strategy of 
introducing a norm and then minimizing the associated cost function 

d>:=||/-^|| 2 = ||/|| 2 -2(/,^) + ||^|| 2 (3) 

can then be used to determine the unknown source coefficients. In order to consider not only fits like ([1]) 
and but other general forms as well, suppose that ip is a linear combination of general basis functions: 

N k 
k=l 
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For the real setting, the values of /j, k that minimize $ = ^(/ifc) can be found simply by solving the linear 
equation set that results from setting the partials of $ with respect to fi k i to zero (and diving by two): 



^2(B kl ,B k )fx k = (B k >,f), (5) 



k=l 

which can be written more compactly as 



T> = A (6) 



where T denotes the matrix whose elements are Ty tk := (B k /, B k ), A the vector whose element are A k ■= 
(B k , /) and /x the vector whose elements are fj, k . 

For the complex setting, when complex conjugation is applied to the first factor in an inner product, 
rather than the second, it is said to be in left conjugate form (LCF) [4] and it is a simple matter to show 
that ([I]) also specifies the solution that minimizes the cost function <J> specified by (J3J) 0]. In the complex 
setting T is obviously Hermitian. 

In the real R n setting for n > 2, it is assumed that / and g are harmonic over f2, which may be an 
unbounded domain, and that the inner products are specified in terms of Dirichlet integrals so that |/| and 
\g\ must fall off to zero sufficiently fast as \X\ — > oo [4J, [5] . Recall that the Dirichlet integral is usually denoted 
by D[/, g] — L. V/ • Vg dV and that this notation was generalized in an obvious fashion in |5j to include a 

weighting function ji = (J,(X) > and to indicate the domain f2 of interest: 

B\f,g,n, 0] := / (iVf-VgdV. (7) 

h 

While only the case \i = 1 is of interest for half-space and D[/, g, 1, f2] = D[/, g], this generalized Dirichlet 
integral notation will still often be used here in order to indicate the domain under consideration. Thus, for 
R™ half-space the primary inner product under consideration will be 

(/, ff ) D/H „ := D[/, .9, 1, FT] . (8) 

Whereas these same approximation considerations apply to any appropriate set of basis functions {Bk} 1 ^^ 
defined over the region of interest f2 ; the interpolation attributes hold only for f2's with certain specific 
geometries of interest and these attributes result from the interaction of the associated norms defined over 
Vt with these particular basis functions. Specifically for an appropriate R™ inner product environment, 
where Q C R™, if the basis functions B k (X) have a point reproducing or replication property, then the 
inner-products of interest involving B k can be written in the form 

(B k , /) = h(P k )f(P k ) , (9) 

where P k G Q and h(P k ) ^ is a fixed function [5]. Corresponding point reproducing or replication 
relationships hold for the complex setting, but for the purposes of discussion just the R™ setting is considered 
for now. [Also although ignored here for simplicity, in general, the possibility of a linear differential operator 
is also considered on the right hand side (RHS) of © — see, for example, (HU).] 

Observe that there are two general situations of interest here. The first is when B k (X) can naturally 
be written as a function of two arguments that are both in fi, say B k = B(X, P k ) where X and P\ G O , 
and the second is when B k {X) can naturally be written as a function of two arguments that are in two 
different domains, say B k = B(X, X' k ) where X £ £1 and X' k £ fl' , which is the compliment of f2. This 
second situation is of special interest in what follows and, further, XL = X' k (P k ) and P k — P k (X' k ) hold 
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for all the particular cases studied here, where P k G [Observe that for reproducing kernel Hilbert spaces 
h(P k ) = 1 and K(X, P k ) ■= B k {X) = B(X, P k ) = K{P k , X), where X and P k G Q and K(X, P k ) is called 
the reproducing kernel.] When a set of basis functions satisfy a point reproducing or replication property 
and T _1 exists then one can show that solutions specified by ([6]) have the interpolation property: 

<p{Pk) = f(Pk) for h = 1, 2, 3, . . . , N k , (10) 

but how can one show in general that |T| ^ 0? 

For the regions of interest, f2, it was shown in [5] that harmonic functions have a strong enough property 
to guarantee that all sets of basis functions satisfying a point reproducing or replication property of the form 
(|9|) are linearly independent. A set of functions with this property is labeled uniformly pointwise independent 
and it can be defined as follows: 

Definition [TJl A set J of scalar valued functions is said to be uniformly pointwise independent over a 
region f2 C K n if for every finite N the following holds: For any bounded set of N arbritrary points Xj G 
and arbritrary bounded constants Cj, there always exists a function / 6 J such that f{Xj) = Cj for 
j = l, 2,3, ...,N. 

As shown in [5] , when a point reproducing or replication property of the form (JS|) holds, this condition is 
sufficient to guarantee that the following two distinct minimum norm conditions are simultaneously satisfied: 

• First, 

II/ — <p\\ is minimized (11) 

directly by ([6]). 

• Second, of all those functions g that satisfy 

g(Pk) = f(P k ) for k = 1, 2, 3, . . . , N k , 
ip, with coefficients fixed by ([6|), has the least norm: 

IMI > IMI ■ (12) 

(Trivially, if a 2 > b 2 then |a| > |6| and conversely, so that ||/ — ip\\ 2 and \\f\\ 2 are also concurrently minimized.) 
For the various DIDACKS half-space settings, since the norms of interest are energy based norms, these two 
minimization criteria imply the associated solutions arc minimum field energy solutions, which has significant 
physical implications in-and-of itself. 

Definition [1] 1 , as well as the associated discussion centering on © through (fT!?)) . extends in a obvious 
way to the complex setting so this extension is automatically presupposed in the sequel. 

Having considered the general framework, consider specific examples of inner products and kernels of 
interest from 4, 5 a , which display a point replication property: 

1. Standard Integral Norm Setting for the Exterior of a Complex Unit Disk: 

Here it is assumed that \z\ > 1 and that / and g have a power series representation of the form: 

= ( 13 ) 

with a.j G C; however, the a/s need not be known explicitly. Then the kernel l/(z — z k ) for < \z k \ < 1 
has the point replicating property [3] 

((z-z k )-\ f) a =p* k f(pl) , (14) 
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where pk — 1/zu and the associated standard integral inner product is defined by 



2tt 

(f>9) a ~±- j [(/(*))*$(*)] de (is) 

6=0 

in the usual polar coordinates (i.e., z = re l& , Zk = rke l6k ). 

2. Dirichlet Integral Norm Setting for the Exterior of a Complex Unit Disk: 

Here it is again assumed that \z\ > 1 > \zk\ > and that, at least in principle, f(z) and g(z) have the 
series form (|13|) . Then it is easy to show that the kernel l/(z — Zk) has the point replicating property 



((z-z k )-\f) n = -(4) /'(Pfe). (16) 



>D 

where again pk := 1/zk and 
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/'(*) := d ~f ■ (17) 
dz 

Here the complex form of the Dirichlet integral inner product is given by 

oo 2?r 

<'■'>»-*//(£)*(£)"*'"• < 18 » 

r=l 9=0 

where, as above, a factor of 1/2tt has been introduced for convenience. (The relationship of this inner 
product to the corresponding R 2 Dirichlet integral D[/, g, 1, f2], including the factor of l/2n, was 
discussed at length in [3].) 

3. Complex Unit Disk Logarithmic Dirichlet Integral Norm Expressions: 

The inner product is the one just studied for the exterior of a complex unit disk and Zk and pk have 
the same meaning. Now, however, one basis function choice to be considered is 

^ fc (s)=ln 2 > (19) 

[Z ~ Zk) 

and for this basis function the point replication condition becomes 

W-fc, f) D = \ f(P*k) ■ (20) 

The basis function choice specified by (|19p is only a special case. In general it is useful to consider the 
paired point formalism introduced in [4] , which has the associated basis function 



Zk(z) := In In T = In k - . (21) 

z - z k z - z k \z— z k J 

Here z' k is called the paired point of Zk and it is an arbritrary complex number, except for the restrictions 
that \z' k \ < 1 and z' k ^ zy for k' = 1, 2, 3, . . . , Nk- [Clearly the choice z' k — for all k results in (fTT)|) .] 
For a general paired-point basis function 

{^f) D = \f{pl)-\f{{p'kY) , (22) 

where (p' k )* is the conjugate involution of the paired point of Zk and is thus specified by taking the 
conjugate of p' k := l/z' k . 
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4. Dirichlet Integral Norm Setting for K 3 Half-space: 

Here © becomes 

k—1 

where vectors in the compliment of the field region are denoted by primed vectors so that £k '■= 
\X — X' k \, where X' k — (x' k , y' k , z' k ) T denotes a typical source point for k = 1, 2, 3, . . . , Nk- Since 
X = (x, y, z) T G H 3 and X' k e H 3 , z > and z' k < 0, so that l/i^ 1 is always bounded for a specified 
X' k . (Of course, it is also assumed that X' k , ^ X' k for all k' ^ k.) The H 3 point replication (or 
generalized collocation) property for the DIDACKS kernel i^ 1 is simply 

D[/,^\l,H 3 ] - 271-/(4, Vk, '4) • (24) 



Higher-order pole and multipole fits 

Other relationships similar to (fl"4|) . (|T6|) . (f20f and (|22ll were given in T|, including corresponding R 2 
point replication conditions. A discussion of analogous relationships for the interior of unit disks was also 
included. Only the complex half-space relationships corresponding to (fTl")) . (flUj) . (|2TJ)) and (j!?2"j) are derived 
here since the other corresponding R 2 relationships can be derived by a relatively straightforward repetition 
of the steps followed in [4j. Further, in the complex setting, as observed in [4], higher order poles can be 
easily incorporated into the fitting function if since 

1 \ (-l) m m! d m ( 1 



dz m \z-z k J (z-z k ) m+1 dz^\z-z k 



(25) 

and when partials with respect to Zk are taken on both sides of (fT4|) analogous closed-form expressions are 
produced for higher order poles. For example, a general fitting function of the form 

N k N m (k) 

k=l 771=1 V ' 

can be easily handled, where the i>k, m's are the different pole strengths of various orders at distinct locations. 
(Here N m (k) is the number of different types of poles at location z k .) Clearly when this form is substituted 
in place of ([T]) into $ := ||/ — ip\\ 2 a and partials with respect to the pole parameters are taken, a set of 
closed-form linear equations results for these various higher-order pole strengths. 

In [S] , point replicating expressions analogous to were given for both the interior and exterior of an R 3 
unit sphere, which in conjunction with the corresponding results for R 2 units disks, leads to the assumption 
that analogous results can easily be derived for the interior and exterior of unit R™ hyper-spheres for n > 3; 
however, it is unclear what the applications of such results might be so they are not considered here. 

With the formalism developed here, just as higher-order pole fits can be easily be done in the complex 
half-space setting, in the R™ half-space setting higher-order multipole fits can be done. First, observe 
that when the partials of @ with respect to the components of X are taken a fitting form results that 
is a combination of higher-order (electrostatic) multipole basis functions. Second, observe that partials of 
fundamental solutions with respect to the components of X and with respect to the components of the source 
coordinates X' k differ by only a sign. For example, in R 3 

!-7- ± ~T~!-7- i ~7- <27 » 

dx k X-X k dx X-X k 
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Third, observe that when = Pk(X' k ) is assumed and the partials of both sides of © are taken with respect 
to the components of source coordinates, then a closed form expression results on the RHS and these source 
partials can be taken inside the inner product on the left hand side (LHS). In conjunction with the two other 
observations, this means that if a DIDACKS point replicating inner product exists for fundamental solutions, 
then closed-form inner-product expressions also exist for all higher-order multipole basis functions, such as 
dipoles and quadrupoles, and that interpolation fits for these higher order multipoles can easily be done. 
(For further details see [4] and [5].) Due to fairly obvious connections between inverse- multiquadric radial 
basis function fits and higher-order multipole fits, the M. n half-space point replication conditions developed 
in Section [5] have clear implications for the reinterpretation of inverse-quadric radial basis function fits. 

Since the obvious connections between inverse-multiquadric radial basis functions and higher-order mul- 
tipoles are invariably ignored (or have generally gone unrecognized) , it is worth explicitly pointing them out 
here since they also imply that the associated radial basis functions fits can also be considered DIDACKS 
fits as well and there are advantages to doing this. These connections are discussed in the next section. 
In Section [31 it is pointed out that a development similar to the DIDACKS half-space one can be carried 
out for the half-space surface norm and this implies that a corresponding inverse-multiquadric radial-basis 
function reinterpretation can be carried out for this norm as well. Thus, from (fTTj) and (fl2j) . this also implies 
that there are associated surface integral minimization conditions that hold implicitly for certain radial basis 
function fits. 



2 Reinterpretation of Radial Basis Function Fits 

This section gives a sketchy overview of techniques for reinterpreting radial basis function fits as DIDACKS 
fits. The same issue is briefly addressed from a different perspective for a different norm setting in the next 
section. 

Consider a typical radial basis function fit to prescribed function values fk ■= f{Qk) for k = 1, 2, 3, ... , Nk , 
where Qk G R™- As normally performed, this R™ fit can be boiled down to the following three steps [2j [6]: 

• First, an appropriate radial basis function shaping function <£(|X|) is chosen, where X £ R™. 

• Second, it is assumed that the candidate fitting form is a linear combination of shifted shaping functions, 
so that 

<p(X) = Yl »k<f>{\X-Qk\) , (28) 
fe=i 

where Hk G R- 

• Third, and finally, the coefficients fik are determined directly by solving the point matching condition 
that corresponds to (TTUj) : 

p(Qk) = fk for k = 1, 2, 3, . . . , N k . (29) 

One commonly used type of shaping functions are the inverse multiquadrics, which have the form {see, for 
example, [2] or equation (5.3) of [6]}: 

«r) = ^hy (30) 

where r = \X\, f3 > n/2 and a is called the shape parameter. To indicate how these inverse-multiquadric fits 
can be reinterpreted as DIDACKS fits, several specific examples will be considered — from these examples 
the general strategies to be used in this reinterpretation process should be obvious. 
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Before proceeding, notice that for both DIDACKS fits and radial basis function fits, that if individual 
basis functions are defined by bBk or b<p for some constant &, rather than Bk or (ft, then the same final overall 
fitting expression results when a fit is carried out. The constant b will be called a basis function scaling 
factor. 

The value (3 = n/2 is also occasionally used in the form (|30|) to define the shaping function of choice. 
When this particular shaping function is used to perform R 2 radial basis function fits, ([28]) becomes 

y 1/a (X):=X: l = (31) 
fe=l ^\X-Q k \ 2 + a 2 

which corresponds to an R 3 half-space point source fit of the form 

(p(x, y, z) = > ; (32) 

provided z' k — constant, Qk := (x' k , y' k ), z := \z' k \, a := 2\z' k \ and <pi/%{x, y) :— <f{x, y, a/2). [Notice that 
here, as well as in the sequel, this reidentification process really only requires that z = constant in {32]).] 
Next observe that a dipole oriented along the z— axis in R 3 yields a basis function of the form 

d 1 z — z' k 

k (x, y, z) = --^============ = 

and thus if a basis function scaling factor of l/(2z' k ) is introduced, the same reidentification process that 
was used to match up (f32|) with (j3Tj) can again be carried out for (3 — 3/2, when n = 2 in (f30|) . Taking 
higher-order partials results in like correspondences for (3 = 5/2, 7/2, /9/2 etc. In all of these cases, observe 
that from (JMJ) , (jlip and (fT2"]) , one can infer that the associated radial basis function fits also satisfy analogous 
energy minimization conditions — although the exact interpretation and forms of these corresponding energy 
minimization conditions is left to the reader. 

Next considered the generalization of this reinterpretation technique from R 2 radial basis function fits 
to R™ radial basis function fits. As a first step, for (3 > 1/2, notice that another obvious that way to 
obtain a correspondence between R 2 radial basis functions and fundamental solutions exists for half-space 
fundamental solution fits in R m , for m > 3. Thus consider the R m form 

F = — — , (34) 

\X-X' k \™~* 

which was introduced in @. Let X' k = (x' k , y' k , 0, 0, • • • , 0, z' k ) T , then using the same reidentification 
technique as before it is clear that for ip = tp(x, y, 0, 0, • • • , 0, |^|) T , the associated F corresponds to an R 2 
inverse multiquadric with (3 — (to/2) — 1. This same reidentification procedure can also be carried out for 
higher dimensional radial basis function fits in an obvious way and thus corresponds to one of the desired 
general reinterpretation techniques. Also, as before various partials of F can be taken, which yields another 
reinterpretation procedure for n > 2. As before, the interpretation of the associated energy minimization 
principles is left to the reader. 

There is one significant difference between the typical reidentification process epitomized by (|32p and 
(f3"3")l : Namely, when a partial derivative is taken, say with respect to z as in (|33p. the resulting DIDACKS 
fit corresponds to matching point data for df /dz and not / itself. This, in turn, means that the associated 
inverse-multiquadric radial basis fit must be considered a DIDACKS fit that minimizes the difference of 
the appropriate multipole form and corresponding integrals of /, which is a significant conceptual issue in 
itself. However, there is a fairly direct way of injecting the fitting form (|33[) into surface integral dual-access 
collocation-kernel fits, as briefly mentioned in the next section. 
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3 MJ 1 Half-space Surface Integral Theory 

For the complex setting, (115[) implies that a standard (i.e., surface) integral norm point replication theory 
exists for C disks, while (|46[) implies that a similar development holds for C half-space. Furthermore, from 
the complex to real norm mapping given in [3], it is easy to obtain corresponding R 2 half-space relationships. 
Section [6] develops an analogous theory for R" half-space (for n > 2). In the complex setting the associated 
functional space setting was labeled a standard integral dual-access collocation- kernel space (SIDACKS). For 
the R™ half-space setting the same acronym will be reused to mean a surface integral dual-access collocation- 
kernel space. It is perhaps surprising that such a theory exists, but it is even more surprising that it 
has been (completely?) overlooked in the literature given its direct and obvious link to Poisson's integral 
formula for R™ half-space. The resulting surface inner products between dipoles that are oriented along 
the upward pointing unit norm direction, such as (|33|). and some given R™ harmonic function / yield a 
result proportional to a point evaluation of / itself and not the partial of /, which implies an alternative 
reinterpretation of inverse-multiquadric radial basis function fits. Moreover, this reinterpretation obviously 
implies that corresponding surface integral norm minimization conditions hold, in complete analogy to the 
energy-norm case. Again this actual reinterpretation of inverse multiquadric radial basis function fits as R" 
SIDACKS fits and the associated minimum norm implications are left to the reader to sort out. 

There is one other significant problem area where SIDACKS theory is of special interest. This is the prob- 
lem area of downward continuation in R™ half-space. There seems to be no existing appropriate physically 
motivated mathematical framework for rigorously handling this type of problem (although various mathe- 
matical formalisms do exist). The general issue of using a (weighted) energy- norm approach was raised in 
Appendix B of [5. and it will be followed up in a subsequent article in this series. Here downward continua- 
tion will be considered from a different perspective for half-space geometries in R™ for n > 2, although it is 
usually only the R 3 half-space case that is of interest. In general, the R n half-space downward continuation 
problem is that of reconstructing the value of / on some surface closer to <9H™ from specified values at a 
surface that is further away from dH n and labeled E. Here it will be assumed that the surfaces in question 
are parallel to <9H™. Furthermore, for convenience, the surface where / is to be reconstructed will be taken 
to be coincident with dH n itself. 

Let ip denote the full reconstruction that is to match / on the surface E. First, observe that since the 
extremes of a harmonic function occur on the boundary, what one would ideally like to minimize is: 

$00 := ll^-ZIll wh ere (35) 
||/||^:= Sup |/(P)| 2 for all PedK n . (36) 

Notice here that 11/H 2 ^ is the square of ||/|| := Sup \f(P)\ for all P G <9H n and that it is somewhat hard 
to deal with computationally, so it is expedient to seek out an alternative minimization form. Reasoning 
along similar lines, it would seem that the second best form to minimize would be 

*„/H» :=\\v-ft /nn where (37) 

(/. Sl/H- : = / fa dX (38) 

and where d X is defined in an obvious way in Section[5] Realistically, in general it is necessary to reconstruct 
/ from data values that are specified over only part of E and where noise is present in the data that is available. 
The issue of noisy data is ignored here, but if appropriate inner products can be evaluated from sufficiently 
accurate available data, then f> can be used to do downward continuation. In Section [6l it is shown that 
two types of surface inner products can be theoretically evaluated. One where the inner product involves 
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dipoles and is naturally associated with point evaluation on E and one where the inner product involves 
fundamental solutions, but that requires information about / that is in the region above E [in particular, it 
requires line integrals of / outward along the (upward pointing) normal direction of <9H™ from some specified 
point of E] . 



4 Point Replication Conditions for Complex Half-space 

In this section the replication conditions for half-space corresponding to (|14|). (fl6|) . (|20|) and (|22|) are derived. 
There are three general strategies for obtaining these corresponding relationships: 

(a) Take the appropriate limits of (JU) , dlK (EOJ) an d (J22) by considering what happens when z k ~ (0, i) 
and then translating to new coordinates. 

(b) Develop entirely new relationships using standard analysis techniques. 

(c) Recast existing RKHS complex half-space (or half-space boundary) kernels into dual-access collocation- 
kernel form and then derive the corresponding point replication conditions. 

The last approach was addressed in [4] where, for example, various specific existing kernels were analyzed 
for both the interior and exterior of complex unit circles and unit disks; moreover, it was also pointed out 
there that a general correspondence exists between a Bergman kernel for any specified complex domain 
and a complex Dirichlet integral reproducing kernel defined over the same domain, so that Bergman kernel 
results for half-space can be recast as Dirichlet integral dual-access collocation-kernel results. Nevertheless, 
unless the study of Bergman kernels is the primary focus, approach (c) is not really as enlightening as cither 
approach (a) or (b), so it will not be followed up here. (As an aside, it is an unstated implicit goal of the 
DIDACKS programme to develop all the associated kernels independently and then recast them in RKHS 
form in order to verify that no significant RKHS possibilities have been overlooked.) This leaves a choice 
between approach (a) and (b). Although both approaches are instructive and there are some subtleties 
involved in carrying out either one; the second approach is ultimately the more satisfying and convincing 
one so it is the main one opted for here. 

Clearly the unmodified series form (fT5|) is inappropriate for the half-space case due to its behavior 
at the origin of H c . Thus for both the standard integral complex half-space setting and the Dirichlet 
integral complex half-space setting, it is assumed that all admissible analytic functions / have a power series 
representation of the form 

oo 

in H c for bounded d € H^ and / that tails off "sufficiently fast" as \z\ — > oo. (Here "sufficiently fast" should 
be obvious in the sequel from the context.) 

Looking ahead, for each point replication condition it is necessary to consider, not only the definition of 
the inner product, but of the fundamental reproducing kernel, Fk, as well. The point is that although it 
may be assumed, for example, that Fk has the form of a simple pole, for convenience it may be necessary to 
include an overall basis function scaling factor so that for simple poles the general kernel form assumed will 
be F k =b/(z- z k ) for z k E H^ . 
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Standard integral norm for the complex half-space setting 



For z G H c the natural analog of (fT5j) to consider is 



(9,1)^:=^ I Mz)Tf(z)} 



y=o 



(40) 



x— — oo 



Here the integral 
that 



/llo-/2 mus t be bounded. Also observe that from the assumed series form (j3T))) it follows 



lim 



[(<?(*))*/(*)] 



9=0 



Next consider the Cauchy integral formula 



/(*) 



c 



(41) 



(42) 



where / is analytic in the simply connected region bounded by the closed curve 7 and £ is in this same 
region. Taking (|4"T]) into account and choosing 7 in (|4"2"|) to be the path from — ii to i? along the x— axis and 
then from 9 — to 8 = tt along a circular arc of radius R, immediately gives 



/(C) = ^ 

2m 



c. 



(43) 



for admissible functions / when R — > 00. The exact choice of the form for i 7 ^ remains, but a comparison of 
(|40| and (|43|) is suggestive. 

In [3] it was pointed out that l/(z— Z}~) corresponds to an R 2 dipole term oriented along the x— axis and 
that ij(z— Zk) corresponds to an R 2 dipole term oriented along the y— axis. From the orientation chosen for 
the half-space H c itself (i.e., y > 0), it is obvious that fundamental kernels should correspond to R 2 dipoles 
oriented along the y— axis. Thus introducing a basis function scaling factor of b = i results in the choice 



F k (z) 



Zk 



Substitution of this expression for g into (|40|) immediately yields 

f(x, 0) 



(i(z-z k )-\ /) 



4) 



d x 



(44) 



(45) 



(As an aside, since the term "fundamental solution" is not normally used in the complex setting, the term 
"fundamental kernel" is used here to mean a simple pole; alternatively, in the M 2 setting logarithmic sources 
are fundamental solutions so from the natural correspondence implied by standard harmonic completion [4] 
one could take the term "fundamental solution" in C to be a logarithmic term, but l(z — z k ) is much more 
ubiquitous than logarithmic terms are in C so a new phrase with a new meaning is introduced here to simply 



by-pass this issue.) Introducing the half-space relationship q k 
desired result: 



(F k , f) 



<t/2 



f{ik) = m) . 



and comparing (145(1 with (|4"3")l gives the 

(46) 



Here q k is the half-space analog of the unit disk interior and exterior conjugate involution point of z k : 
p k = which should be obvious from a consideration of the limiting process mentioned at the first of 

this section in connection with approach (a). 
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Dirichlet integral norm for the complex half-space setting 



The natural half-space analog of (fl8|) to consider is 



DC 



X— — 00 y=0 

First, consider the choice of the overall basis function scaling factor. In (|4"4")l the choice b = i was justified 
on the grounds that it corresponds to a M 2 dipole oriented along the y— axis. This choice could have been 
argued in perhaps a more precise way by looking at the limiting process of (fl4)) mentioned at the first of this 
section in connection with an alternative derivation process [approach (a)]. In this limit for (|14|) ~ i, so 
it is obvious that an extra factor of i enters in the final limit from the RHS of (114|) , which corresponds to a 
choice of b = i. When the same limiting process of (|16p is considered, then (pi) 2 ~ —1 occurs in the limit of 
the RHS of l|16p and thus, aside from a possible overall sign difference that can be accounted for later, there 
is no need to introduce a basis function scaling factor. Hence for convenience the choice 

«W := ^ (48) 
is made here and the inner product of interest becomes 

oo oo 

(F fe , /) :=-! / / - f{z) dydx, (49) 

u l 2 2vr J J (x-iy-qk) 1 

x— — oo y=0 

where, as before, q k — zj* for the half-space setting here [see the comment after (|46p]. 

Observe that due to the factor of — iy rather than iy occuring in the denominator of the integrand on 
the RHS of (I49|) . the Cauchy integral formula cannot be applied directly, and, furthermore, it is hard to 
transform this factor of — iy away while still preserving the analytic character of the rest of this integrand; 
hence, a certain amount of formal manipulation seems to be needed. Thus, first let 

oc 

K{x)--= I 7 ! ' {Z) s 2 dy (50) 

so that (j49|) becomes 

^f) D /2 = ^ J K{x)dx. (51) 

X— — OO 

Since 

d r _ % , if 



dy(x-iy-q k ) (x-iy-q k ) (x-iy-qk) 2 
and from the Cauchy-Riemann equations 

dy dx 
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(i.e., if g = u + iv then [du/dy] + i[dv/dy] — — [dv/dx] + i[du/dx]) it follows that 

df 



K{x) 



/'(*) 



(x-iy- q k ) 



<).!■ 



d_ 

dy 

v=o 

if'jx, 0) 

(x-qk) J {x-iy-qk) 
v=o 



dy + i 



y=0 



(x-iy- q k ) 



dy 



dy 



(52) 



Likewise 



dx 



r 



so that 



JC{x) = - 



o r = 

dx(x-iy-q k ) (x - iy - q k ) (x - iy - q k ) 2 

CO 



d_ 

dx 



{x-iy- q h ) 



dy- 



dx 



(x-iy-qk) 



dy 



(53) 



y=0 y=0 
Adding to the extreme LHS and RHS of (|5^)) and dividing by two yields: 



JC{x) 



if'(x, 0) 1 f d 



2(x — qk) 2 J dx 

v=o 



f'(z) 



{x-iy-qk) 



dy. 



Thus after the indicated evaluations of the resulting double integral are performed the following desired 
expression results 



fC(x) dx 



f'{x, 0) 

dx = 

(x - qw) 2ni 



f'(z) 
z - q k 



dx = -nf'{q k ) 



(54) 



y=o 



where the Cauchy integral formula given by (|43| was used the evaluate the RHS here. 
Substituting this result into (fSTj) immediately gives the end result: 



(( z - Zk )-\f) D/2 = lf'{ qk) = l±- 



(55) 



Complex half-space logarithmic Dirichlet-integral inner-product expressions 

The inner product is the same H c one just analyzed; furthermore, z k and q k have the same meaning so 
that q k — zl with lm{q k } > 0. Now, however, the following basis function will be considered, which is the 
obvious analog of (fT5|) : 

^( 2 ):=ln-^-. (56) 
[z - z k ) 

The following analog of (|2"Tj) will also be considered: 

&(*) := In - In = In ( Z —^) , (57) 

z - z k z - z, \z — ZkJ 



13 



which is a generalization of (|56p . Here Im{z^,} < and it is useful to introduce q' k :— (z' k )* . Although the 
factor of i occuring on the RHS of (|56p looks strange and it occurs because ipk(z) is not translation invariant 
[while £k(z) is], it is expedient to derive expressions for it first. 
Observe that 

dipk 1 1 ,-„s 

—, — = ; ~ - 7 7 • (58) 

dz (z + i) (z — Zk) 

Thus from (|47p . the inner product expression of interest is 



oo oo 

1 



X— — OQ y—0 



y-i x-vy-qu 



f'(z)dydx. (59) 



For the sake of variety, consider a Wirtinger calculus based derivation approach [4]. First, since the 
Wirtinger calculus is based on the concept of conjugate variables, recall that the notation (z, z) was reserved 
for this pair of variables in [4] so that, in general, / = f(z, z). Since / is analytic it does not depend on z 
and 

dz d z 

so that the usual Wirtinger derivatives {c.f., equation (1) of [4]} 

d 1 f d d 
— i- 



dz 2 \dx dy 
can be used to reexpress ([59^1 : 



and 



— i- 



d x (x — iy — a) (x — iy — a) 2 (x — iy — a) dx 
.d f I . 1 Of 



dy (x — iy — a) (x — iy — a) 2 [x — iy — a)dy 

oo oo 

^ f) = J_ f f 1 (9 .d \( f / \ dydx 

D / 2 2n J J 2\dx dyj\x-iy-i x — iy—qh) 



X — — OC y—0 



Obviously 



and 



d \f f I 



d x J \x — iy — i x — iy — qu 



dx = 



u=0 



r-- — ) dy=-rJL__L_ )/(X)0 ) 

dyj\x-iy-i x-ty-qkj \x - 1 x - q k / 



(60) 



^^W = i / / {■ X -ly— l -x- i 1 y-q k ) K^-^) dVdX - (61) 

X — — OC y — 

Now for any constant a E C 

d f -f 1 df 



(62) 
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( ^' / W = ^ y 2\irT-T^) dx > (63) 

a;— — oo 

which involves the same form of integrals that occured on the RHS of (|^5j) or ([33]). Hence, from ([^5)1 

(V-fc, /) D/2 = | /(ff*) - | /(i) • (64) 
Retracing the steps in the above derivation it is obvious that 

(Ck,f) D/2 = \f(qk)-\f(q' k ), (65) 
in accord with the complex half-space limit of (f2"2"|) . 

5 Dirichlet Integral Replication Expressions for W 1 Half-space 

This section derives R" half-space results analogous to ([24]) . As previously indicated, the R 2 half-space 
inner-product expressions can be obtained in a fairly straightforward manner from the Dirichlet integral 
complex half-space results, so only the case n > 2 will be considered. 

First, consider the question of overall notational conventions. While there are several older reference 
works that describe general R™ potential theory, such as [3J, it is useful to stick closely to more modern 
treatments since several notational issues arise in R n that do not arise in R 3 and it is a good idea to have 
a reference that is readily available. Since [JJ gives a fairly broad treatment of R™ potential theory at the 
introductory level, since it contains all of the necessary results, since it is in print and widely available and, 
finally, since it is well written it will be used as the primary reference here; moreover, to circumvent the 
issues just noted, the notation adopted here will be the same as that used in PQ except for a few indicated 
or obvious exceptions. 

Next, consider the task of reexpressing the Dirichlet integral given on the RHS of (|8|). Three mathematical 
tools are useful in this endeavor: 

• The divergence theorem {[TJ p. 4]: Equation 1.2} 



V -W dV = Jn-WdS (66) 
n an 

where W is a vector field in R™ [W = W(X) for W and X £ R ra ] and n — n(X) is the outward pointing 
unit normal vector associated with the boundary surface <9f2, which has a differential hyper-surface 
area of dS. 

The product rule for the Laplacian given by {[TJ p. 13]: Equation 1.19} 

V 2 (uv) = uV 2 v + 2Vu- Vv + vV 2 u . (67) 

Green's identity, which can be written as {[TJ p. 4]: Equation 1.1} 

(uff - uff ) dS = J (u\7 2 v - v\7 2 u) dV . (68) 
on n 
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If V 2 </> = V 2 -0 = holds over all of fi, then from (|6T)) 

and it follows from (|66[) that 

D[0, V, 1, 0] = i /" V • V(#) dV=^Jn- V(#) dS=\J ^(#) . (69) 
a an an 

Further it follows from (|68[) that 

f&*S = f*%*S (70) 
ao ao 

so that finally 

d[^, ^ i, n] = /" vH ■ (7i) 

Consider the specialization of (fTTj) to H™. Following [I], let Z G H n C K" denote the vector variable of 
interest and reexpress it as a combination of two parts: For n — 1 > j > 1 let a^- = and for j — n let 
y = z„, so that Z = (zi, z 2 , z 3 , ■ ■ ■ , z n _i, = (art, x 2 , £3 • • • , x„_i, y) T = [X T | y } T where X e R n_1 
and X := (x\, X2, £3 ■ ■ ■ , x„_i) T . {Here Z was written as [^ T | y] T rather than in the equivalent form 
(X T , y) T in order to emphasize the fact that a block matrix partition is implied.} Further, let dX — 
dx\ dx2 dx3 ■ ■ ■ dx n -\ and let 



00 oc 

/ ■ ■ / 



Thus, since d/dh = —d/dy for the half-space H n given by y > 0, it follows from ([71]) that 



D[</>, V, 1, H n ] = - J ^ dX . 



(72) 



The strategy will be to compare the results from an evaluation of (|72| for a suitable choice of <p with 
Poisson's integral formula. Temporally leaving aside the choice of (j> to be made, consider Poisson's integral 
formula for 1" half-space. From [TJ p. 9]: 

1>W= J Pn^(t,0)dt (73) 

where t G R™ -1 is restricted to <9H™, just as X £ ]R™ -1 was and dt is defined in complete analogy to dX. In 
(|73J) P H = P H {Z, t ) has the form [TJ p. 145] 

^ (2, t* ) = ^ (74) 

(|Jjf-t |2 + y2)»/2 

where [TJ p. 145] 

C " = ^5) (75) 
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with [H p. 245] 



V{B) 



T n/2 



r(| + i) 



so 



ip(z) = c n y 



(\X -t\ 2 +y 2 ) n / 2 



dt . 



(76) 



(77) 



Next consider the choice of basis function and associated global basis scaling factor. From the discussion 
surrounding ([2]) it is clear that the general form of interest is 



F(Z,S) = 



d„ 



Z-S\ 



(78) 



for some appropriate basis scaling factor d ni where Z £ H n and S £ H™. For example, with this choice of 
fundamental solution basis function, @ directly becomes 



<P = 



N k N k 

V = » k F(X,X' k 

^ Y - Y n-2 J ^ y k 



(79) 



k=l 



in the notation used there. Although Q] and [3] assume d n < 0, to simplify matters here it is useful to take 
d n > 0. The negative of the basis function scaling factor used in [1] p. 193] leads to the choice 



d n = 



(n-2)nV(B) 2(n-2) ' 



(80) 



which is the one adopted here. 

In analogy to the block matrix partition of Z into X and y, let S = [t T \w] T where w € K is negative. 
Then by a straightforward evaluation 



0_ 

Oy 



F(Z,S) 



Cn (V - W) 

' 2 \Z-S\ n 



(81) 



Substituting 4> = F(Z, S) into ((T2j) thus yields 



D[0, </>, 1, FT] = ^ 



(y-w)j> 
\Z-S\ n 



y=0 



dX = ^\w\ 
2 1 1 



{\X - i\ 2 + w 2 ) n /< 



dX 



(82) 



Introducing P := [t T \ — w 
obvious change of variables) 



\ T so that P e H" and comparing (|82]1 with (|TT[) immediately gives (after an 



D[F,V, 1,H"] = -^(P). 



(83) 



For n = 3, V(B) = 4tt/3 and d n = d 3 = 1/(4tt) so that F(X, S) = 1/(4tt|X - S\) and 1/4 = 
4<kF(X, XjQ. Thus JHH) gives 



D[^\ V, 1, H J ] = 2^(P fc 



(84) 



in accord with 
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6 Surface Integral Replication Expressions for M. 11 Half-space 



As discussed in Section [3] this section derives replication expressions for the 1™ half-space surface integral 
setting (i.e., SIDACKS). Many of the results and much of the notation from the preceeding section are reused 
here. Thus using the notation introduced prior to (|71[) . the inner product of interest here is 



(/, 9) ara n ■= J [fg] 



dX 



y=o 



(85) 



in accord with l|38jl . 
From (1701) 



g^dx 



and from ([72)1 

Which from l|83p implies that 



(<?, df/dy) = -D[g, f, 1, IT] 



(F, df/dy) = --f(P) , 



(86) 



(87) 



where F is given by ([75]) . From ([B"6")l this can be rewritten as 

(F, df/dy) = (dF/dy, f) = (F y , f) 



:f(P) 



where 



F„ 



OF 
dy 



(89) 



(90) 



The last two expressions in (|89)) can be used to perform a SIDACKS dipole fit directly and, as such, can be 
used to reinterpret inverse-quadric radial basis function fits as previously noted; moreover, these expressions 
can also clearly be used to perform downward continuation. 

In many cases it is also possible to do SIDACKS fundamental solution based inner product fits. From 



(91) 



(dF/dy, /) ct/h „ = -(dF/dw, /) CT/Hn =--^(F, f) a/nn = -\f(P) 
Let p y denote the last component of P. Since p y — —w, (|9ip can be rewritten as 



(92) 



which, after a change of dummy independent variable from p y to p' y , can be integrated from p y to oo to 
obtain 



1 



f(P) dp' 



(93) 



P' V =PV 



This last expression can be used to analyze downward continuation, either analytically or numerically. 
Consider the 1R 3 half-space case. First observe that (£/., ^fc') CT m 3 can ^ e eva l ua ted in closed form from the 
RHS of ([93]) so that the associated T matrix can be easily evaluated. Second, for specified f's the inner 
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product (/, ^fc) CT ^ H 3 can also generally be evaluated either in closed form, in terms of a series expansion or 

numerically. For example, in general if a standard M 3 harmonic Fourier integral representation for f(X) is 
substituted into the RHS of (f9"5|) and the indicated integral is performed then the result seems to imply some 
sort of rule-of-thumb, which adjusts R 3 downward continuation effects by a wavelength factor. Alternatively, 
a general numerical analysis of R 3 downward continuation algorithms could be undertaken by using 
for various given mathematical representations of / since a numerical evaluation of the RHS of (|9li|) can 
obviously be performed in this case. An analysis here might consist of comparing these types of results 
with results from a standard downward continuation formalism and with results based directly on <j>oo as 
specified by (|35j) . The last point is that for a given / and specified set of basis functions, the actual (source) 
coefficients that yield a minimum of $00 can be found by employing non-linear least squares algorithms. 

It is also perhaps worth noting here that a if a standard Fourier integral harmonic expansion (for E's of 
infinite extent) or Fourier series harmonic representation (for E's of finite extent) is assumed then fits based 
on (|93p can be performed for noisy incomplete data specified on E (usually either some sort of point or track 
measurements). For example, suppose that a collection of noisy point measurements of f(x, y, z) is available 
for x G {—L, L), y G (— L, L) and z — z > and a reconstruction of f(x, y, 0) [i.e. <p(x, y, 0) w f(x, y, 0)] 
is desired where x G (-L 1 , L'), y G {—L 1 , L') with L' < L. Let the associated harmonic Fourier series be 
written symbolically as (the complete form of the harmonic series can be found in most texts that cover 
electrostatics): 

3^(x, y, z) = {a(k x , k y ) cos(k x x) cos(k y y) + other sin and cos terms} e V k - +k l , (94) 
00 

- 1 / 7(x, ,z)d, = iV M*»,fc») cos(fc x x) cos(M + ■■■} (95) 

2 J 2 . /fc2 I fc2 



SO 

00 

1 f 1 ^ {a(k x , k y ) cos( 

~ n / 2{X, y, Z) d,Z = - > j y - ; 

z=z Y V 

Downward continuation can then be carried out by, for example, following the following three step process: 

(A) Determine the unadjusted expansion coefficients [a(k x , k y ) etc.] by fitting them to the noisy data itself 
by doing a linear least squares (LLSQ) fit to the data. 

(B) Use a frequency domain based procedure to adjust of these coefficients for the presence of noise, so 
as to obtain a nose compensated surface reconstruction of / for E. Let a(k x , k y ) denote the resulting 
best unbiased estimates of the new coefficients. 

(C) Do a fundamental solution downward continuation that is based on inner product evaluations given by 
the RHS of (|95|) . where a(k x , k y ) is to be substituted for a(k x , k y ). 

Several comments are relevant here. First steps (A) and (B) were successfully implemented by the author 
in the early 1980's as part of the analysis for an airborne gravity gradiometer survey system project. (The 
end goal of the project was to obtain surface gravity estimates.) While the details of the implementation do 
not concern us here, it is perhaps worth noting that to successfully implement step (A) requires that certain 
inherent pitfalls be recognized and overcome, but once these underlying issues are recognized the solutions to 
them arc fairly easy to come up with. {The main one arises from the interaction of the Gibbs phenomenon 
with the implied checker-board like periodic repetitions l J{x ± 2L, y ± 2L, z) = ${x, y, z), which produces 
unwanted overtones (consider, for example a one dimension Fourier series expansion of a saw-tooth pattern). 
Not only must long term biases be removed, but the specified data should be edge tapered in order to remove 
all of these unwanted edge induced frequency effects.} Part of step (C) requires that an appropriate grid and 
depth for source locations be selected. It is at this step that some understanding of the interplay between 
source positions, "source regularization" and downward continuation enters. (In Appendix B of [5] it was 
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indicated that an understanding of downward continuation is, in some sense, tied to what is known, or is 
reasonable to assume, in the way of source information or statistics. (Various types of "source regularization" 
will be discussed in a subsequent article dealing with inverse source estimation implementation issues.) 
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